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Abstract 

In this paper, we present our study on the rj photoproduction based 
on the chiral quark model. We find that quark model provides a very 
good description of the ry production with much less parameters, and the 
threshold region is not a reliable source to determine the r/NN coupling 
constant due to its strong dependence on the properties of the resonance 
Sii(1535). We suggest that the systematic data in Ei a b = 1.2 ~ 1.4 
GeV region may help us to determine the rjNN coupling constant more 
precisely. The structure of the resonance 5n(1535) is discussed, we find 
that the recent data from Mainz group bring the helicity amplitude 
much closer to the quark model prediction. However, more studies 
need to be done to understand the large rjN branching ratio of the 
resonance 5n(1535). Our results show that the quark model is a very 
good approach to study the underlying structure of baryon resonances 
from the meson photoproduction data. 



PACS numbers: 13.75. Gx, 13.40.Hq, 13.60.Le,12.40.Aa 



1. Introduction 



In our previous publication [|]], a framework based on the chiral quark model 
to study the meson photoproductions is developed. It starts from the low 
energy QCD Lagrangian|| so that the meson-quark interaction is chiral in- 
variant, and the low energy theorem in the threshold pion-photoproduction[]3| 
is automatically recovered with a proper treatment of the center of mass 
motion 0. By treating the pseudoscalar mesons as Goldstone bosons that in- 
teract directly with quarks inside baryons, the quark model provides an unified 
formalism for all s- and u-channel resonances, and the number of parameters 
used in the model are dramatically reduced. In principle, only one parameter 
is needed for all resonances that contribute to the meson productions. This 
marks significant advance from the traditional theory, in which the effective 
couplings among the hadrons are used so that each resonance requires one 
additional parameter. It also makes it possible to provide a consistent calcula- 
tion of the meson-photoproduction beyond the threshold region. Perhaps more 
important, it provides an unified description for the pseudoscalar meson pho- 
toproductions and highlights the dynamic role by baryon resonances in each 
processes. In this paper, we extend our investigation to the rj photoproduction. 

There are many features of r\ mesons that make the rj photoproductions 
unique in the quark model. The r\ meson is an isospin zero state, thus only the 
resonances with isospin 1/2 contribute in the s and u channel. It is also a charge 
neutral particle so that the contact (seagull) term[Q] that plays a dominant role 
in the charge meson production does not contribute, thus enhances the role 
of resonances. Moreover, because the mass of the resonance Sii(1535) is just 
above the r\ N threshold where the S-wave is dominant, the rj photoproduction 
in the threshold region provides us very important probe to the structure of 
the resonance 5*11(1535). Thus, there have been considerable theoretical and 
experimental interests in studying the rj photoproduction. New experimental 
data for the rj photoproduction in the threshold region from Bates0, ELSA[0], 
and Mainz p| have been published recently. In particular, the data from the 
Mainz group provide more systematic behaviour of the r] production in the 
threshold region, which has better energy and angular resolutions, thus enable 
us to study the properties of the resonance S , n(1535) more precisely. Therefore, 
it is very interesting to note that the helicity amplitude A P i extracted from the 
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new Mainz dataH] is much closer to the prediction of the quark model P, 



T0|| . On theoretical side, the theoretical studies of the r\ photoproduction 



were mostly in the framework of Breit-Wigner parametrizations||ll|| or coupled 



channel isobar models 113]. The recent investigation by the RPI group |T3| 
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has made significant progress in this field, in which the effective Lagrangian 
approach is used so that the properties of the resonance S , n(1535) extracted 
from the data are more model independent, and the number of parameters is 
reduced considerably. 

Traditionally, investigations of meson photoproductions in the framework 
of the quark model have concentrated on the transition amplitudes, in partic- 
ular the helicity amplitudes for the electromagnetic transitions and the partial 
wave amplitudes for the mesonic decays of baryon resonances. These ampli- 
tudes were extracted from meson phot opro duct ion data by the phenomeno- 
logical models, thus less model independent. Instead of relying on the transi- 
tion amplitudes from the phenomenological models, the quark model approach 
enable us to study the structure of baryon resonances directly from the pho- 
toproduction data. Thus, the connection between meson photoproductions 
and more fundamental theories based on QCD can be established. Because 
the transition amplitudes in the quark model have very different energy and 
momentum dependences from those in traditional models, it is by no means 
trivial if meson photoproductions can be successfully described by the quark 
model. This requires that the transition amplitudes in the model have cor- 
rect off-shell behaviour, which they are usually evaluated on-shell. Our early 
investigation [|l| in the Kaon photoproductions has shown that the quark model 
approach presents a much better framework to understand the reaction mech- 
anism of meson photoproductions than many traditional hadronic models, and 
we shall show that the results in the 77 photoproduction are equally encourag- 
ing. 

The paper is organized as follows. The general formalism in the quark 
model for the rj photoproduction is presented in the section 2. We have car- 
ried out three different calculations in the section 3; the first assumes the 
SU(6) ®0(3) symmetry for the baryon wavefunctions so that only one param- 
eter is required to fit the experimental data, the second includes the possible 
configuration mixing effects for the resonances Sii(1535) and Sii(1650) with 
three additional parameters, and the third calculation is concentrated on fit- 
ting the recent Mainz data to extract properties of the resonance S , n(1535). 
Because the data from Mainz group are significantly different from the rest, 
we fit them separately. In section 4, we discuss the structure of the resonance 
5ii(1535) from the 1] photoproduction data in the threshold region, and high- 
light the problems yet to be resolved. Finally, the conclusion is given in the 
section 5. 
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2. General Formalism 



There are two major components in addition to the calculation of the electro- 
magnetic and strong transitions of the baryon resonances in the quark model 
approach, which has been shown to be crucial in deriving the model indepen- 
dent low energy theorem in the threshold pion-photoproductionQ]. First, one 
has to combine the phenomenological quark model with the chiral symmetry, 
this is being achieved by the introduction of the chiral QCD Lagrangianp| 
so that the meson transition operators are chiral invariant. Second, since a 
baryon is being treated as a three quark system, the separation of the center 
of mass motion from the internal motion is important to recover the low en- 
ergy theorem in the threshold pion photoproduction, this has been discussed 
in detail in the Compton scattering 'jN — > 7iV[[5|. 

The differential cross section in the center of mass frame is 

da c - m - _ aea^E 1 + M N )(E f + M N ) |q 



dfl \QsM 2 N 



-\M 



fi\ 
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where a v is the t]NN coupling constant, a e is the electromagnetic coupling, 
^/s = E % + u-y = E* + u v is the total energy in the cm. frame. Generally, 
the 77 transitions between the resonances and the nucleon can be expressed in 
terms of the a v , and no additional parameter for each resonance is required. 
Therefore, the coupling constant has been removed from the matrix element 
Aifi so that it becomes dimensionless. The coupling constant a v is treated 
as a free parameter because of the theoretical issues, such as the Ul anomaly, 
and 77-77' mixing. 

One can write the matrix element Ai fi in terms of the CGLN|§ amplitudes; 

M fi = 3-e (2) 

where e is the polarization vector, and the current J is written as 

t r . ■ , (o- ■ q)(k X a) cr ■ k ct ■ q 

J = f l(T + if 2 _ + / 3 ^^ q + / 4 ^^ q 3 

q k q k q z 



in the center mass frame. The differential cross section in terms of the CGLN 
amplitude is[fT4| 



\M fl \ 2 = Re {|/i| 2 + |/ 2 | 2 - 2cos(6>)/ 2 /* 
sin 2 (^) 



+ 



\h\ 2 + IM 2 + 2fj* + 2m; + 2cos(e)fj; , (4) 
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where 9 is the angle between the incoming photon momentum k and outgoing 77 
momentum q in the center of mass frame. The various polarization observables 
can also be expressed in terms of CGLN amplitudes, which can be found in 



Ref. 14 



The electromagnetic coupling in the nonrelativistic limit is| 



h e = ej 



Tj ■ e 



1 - 



2m n 



[e x k) 



(5) 



and it has been shown Q that the operator h e in Eq. |5] is sufficient to repro- 
duce the low energy theorem for the threshold pion-photoproductions||. The 
corresponding rj transition operator is a pseudovector coupling; 



Tjnr 



A+ 2^ 



Pj 



(6) 



where A corresponds to the center of mass motion, and depends on the mo- 
menta of the initial and final states, and pj is the internal momentum for a 
three quark system. 

Because the 77 meson is a charge neutral particle, the Seagull term that 
plays an important role in the charge meson production does not contribute. 
Thus, the leading Born term would be the nucleon pole term in the S and U 
channels; 



Mi 



-\-ie 6c 2 



10 r, 



+ 



1 



Ef + M N E i + M 



1 



1 



N 



1 



1 



4P i • k 



fJL N \cr ■ e 



2P i ■ k 



a ■ qcr ■ (e x k) (7) 



.Ef + M N E i + Mn , 

where P l ■ k = u^(E l + tf 7 ), /xjv is the magnetic moments of the nucleon, and 
a 2 is the constant from the harmonic oscillator wavefunctions. The matrix 
element for the U-channel nucleon exchange term is 



Mr 



—e 



0O r . 



2Pf ■ k 
1 



LU v k 2 

+ 



1 



Ef + M N 
1 



2 \Ef + M N E l + M 



IN 



+ 



1 



N 



<j ■ e + 



a ■ 



(e x k)<r • qj, 



(8) 



where ■ k = u^i^E^ + |q| cos^). 

The contributions from the t-channel exchange are not included in this 
approach. This has been discussed in some detail in the literature |15[]; if a 
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complete set of resonances is introduced in the s and u channels, the inclusion 
of the t-channel exchange might lead to a double counting problem. This 
may turn out to be an advantage of the quark model approach, since less free 
parameters are needed to fit the data. 

The first excited resonance that contributes to the r\ productions is the 
Roper resonance, Pn(1440). In the SU(6) quark model, its U-channel contri- 
bution is 



M 



Pn(1440) 

1 



-M Pll (i44 )kV 



(Pf ■ k + 5M| li(1440) /2)216m ( 



Ef + M 



N 



a ■ e 



io r , 



Ef + M 



+ 1 



N 



a 2 



uj, n k 2 c[ 2 



a- lE i + M N 
(T ■ (e x k)cr ■ q 



(9) 



where SM 2 ii(im 
given later. 



M 



p n (i440) — A^f^r, and its S-channel contribution will be 



For the excited resonance with higher energy, such as the P-wave baryons, 
we could treat them as degenerate, since their contributions in the U-channel 
are much less sensitive to the detail structure of their masses than those in the 
S-channel. Therefore, we can write their U-channel contributions in a compact 
form; 



Mi 



(Ml + M 



(10) 



The first term in Eq. |10| represents the process in which the incoming photon 
and outgoing 77 meson are absorbed and emitted by the same quarks, it is 

, k ■ q 



+ - 



1 



UJjjUJy 



2m, 



[zA • (e x k) 



1 + 



2m„ 



cr ■ e + 



Q- 



where 



A 



-(Jr. 



+ 



9a 2 m„ 



1 



(Ax(exk))]F(^,P^) 
F^,pf.k + SM*) 
k ' q ,pf -k + 25M 2 ). 



Ae q 
■ ke • qF( 



3a 2 



Ef + M 



+ 1 q- 



N 



'ID 



(12) 



M l + M N Ei+M N , 

The function F(x, y) in Eq. [TT] corresponds to the product of the spatial 
integral and the propagator for the excited states, it can be written as 

M„ 



^ n\(y + nSM 2 ) 



x 



(13) 
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where nSM 2 = (M 2 —M 2 )/2 represents the mass difference between the ground 
state and excited states with the major quantum number n in the harmonic 
oscillator basis, which will be chosen as the average mass differences between 
the ground state and the negative parity baryons so that 5M 2 pa 0.74 GeV 2 . 
The first term in Eq. [TT] corresponds to the correlation between the magnetic 
transition and the cm. motion of the r\ transition operator, it contributes 
to the leading Born terms in the U-channel. The second term in Eq. [11] is 
the correlations among the internal and cm. motions of the photon and rj 
transition operators, this term only contributes to the transitions between the 
ground and n > 1 excited states in the harmonic oscillator basis. The third 
term in Eq. |H] corresponds to the correlation of the internal motions between 
the photon and r] transition operators, which only contributes to the transition 
between the ground and n > 2 excited states. The second term A4jj in Eq. 



10] represents the process in which the incoming photon and outgoing r\ are 



absorbed and emitted by different quarks, and we found that 

Ml = 0. (14) 

This is a direct consequence of the isospin couplings. Eq. [11] can be summed 
up to any quantum number n, however, the excited states with large quantum 
number n become less significant for the U-channel resonance contributions. 
Thus, we only include the excited states with n < 2, which is the minimum 



number required for the contribution from every term in Eq. 11 



For the S-channel resonance processes, the operator A in Eq. | should be 



a --{^WTm^ + 1 ) ci (15) 

in the cm. frame. The calculation of the S-channel resonance contributions is 
similar to that of the U-channel resonance contributions. However, since the 
operator A is only proportional to the final state momentum q, the partial 
wave analysis can be easily carried out for the S-channel resonances. 

In general, one can write the S-channel resonance amplitudes as 

2M R k 2 +q 2 

S — M R 

where y/s = E l + u; 7 = + u v is the total energy of the system, and Or 
is determined by the structure of each resonance. Eq. [16| shows that there 

should be a form factor, e e^ 2 in the harmonic oscillator basis, even in the 
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real photon limit. If the mass of a resonance is above the threshold, the mass 
M R in Eq. [T^ should be changed to 



Mr — » M R (M R — iT(q}). (17) 

T(q) in Eq. |Tj] is the total width of the resonance, and a function of the final 
state momentum q. For a resonance decay to a two body final state with 
orbital angular momentum I, the decay width T(q) can be written as 



r(q) 




with 



and 



|q?l 



^ = r -M^Mf ? _ Mf m 

where Xi is the branching ratio of the resonance decaying into a meson with 
mass Mi and a nucleon, and is the total decay width of the S-channel 
resonance with the mass M R . The function -Dj(q) in Eq. [18| is called fission 



barrier [0, and wavefunction dependent; here we use 

3a 



D l (q)=exp[-^\, (21) 



which is independent of I. A similar formula used in 1=1 tttt and and p-wave 
1 = 1/2 K7r scattering was found in excellent agreement with data in the p 
and K* meson region [0. Generally, the resonance decays are dominated by 
the pion channels, except the resonance 5*11(1535) whose branching ratio of 
rjN channel is around 50 percent. Therefore, we simply set x n = x v = 0.5 for 
the resonance S , n(1535), while x n = 1.0 for the rest of the resonances as a first 
order approximation. 

The operator Or in Eq. [16| can be generally written as 

Or = A [f*tr ■ e + if*(<T ■ q)tr ■ (k X e) + ft* ■ ke • q + fftr • qe • q] (22) 

for the pseudoscalar meson photoproduction, where A is the meson decay 
amplitude and //* (i = 1 ... 4) is the photon transition amplitude. The meson 
decay amplitude A is determined by the spatial wavefunction of resonances and 
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the relative angular momentum between the final decay products. In Table 1, 
we present the amplitude A in the simple harmonic oscillator basis, in which 
the amplitude A depends on the total excitation n and and the orbital angular 
momentum L. The relative angular momentum of the final decay products is 
expressed in terms of the partial wave language in Table 1, in which the 5*, P, 
D and F waves denote the relative angular momentum 0, 1, 2 and 3 between 
the final decay products. The decay amplitude A in Table 1 is the same as 
the expression in Table 1 in Ref. |I(J with g — |/i = and h — Note 

that A has a negative sign, this is consistent with the fitted value for g — \h 



and h in Ref. [10]. 



The photon transition amplitudes f'P in Eq. |22] are written in terms of the 
CGLN amplitudes, which are shown in Table 2. They are usually expressed in 
terms of helicity amplitudes, Au2 and A 3 / 2 , and the connection between the 
two representations can be established. A very important example is the van- 
ishing helicity amplitudes for the transitions between the resonances belonging 
to (70, 4 N) representation and protons due to the Moorhouse selection rule[]l8 



if one uses the nonrelativistic transition operator in Eq. consequently, the 
CGLN amplitudes for these resonances are zero as well. There are 3 important 
negative parity baryons that belong to (70, 4 N) multiplet in the naive quark 
model; 5*11(1650), _Di3(1700) and Dis(1675). However, it has been shown in 
the potential quark model calculation [ITS] that the two states, 70iV( 2 P J 



2o \ i 



M , 



2 



and 70iV( 4 Pjv/) \ > are strongly mixed. Therefore, the contribution from the 
resonance 5*11(1650) to the r\ photoproduction will be studied by fitting to 
the experimental data. Indeed, this will provide us a direct insights into the 
configuration mixing in the potential quark model. 

The CGLN amplitudes for the resonances with total spin 1/2 can be easily 
related to the helicity amplitude A\/2 that has been frequently calculated in 
the quark model. Only the CGLN amplitude f± is nonzero for the resonance 



5n(1535), and this corresponds to Eq multipole transition ||14||. Moreover, 
the amplitudes for the S-wave resonances have the same structure as the 
corresponding helicity amplitude A1/2 in Ref. PB" |, in which the same nonrel- 



ativistic transition operator is used. For the P wave resonances, such as the 
resonances Pn(1440) and Pn(1710), only the CGLN amplitude f£ is present, 
which gives a Mf transition. Notice that the resonances with isospin 3/2 do 
not contribute to the rj photoproduction due to the isospin coupling between 
the r) meson and the nucleon. These results provide an important consistency 
check for the CGLN amplitudes in Table 2. 

If one intends to calculate the reaction beyond 2 GeV in the center of mass 
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frame, the higher resonances with quantum number n = 3 and n = 4 must be 
included. Instead, we adopt an approach that treats the resonances for n > 3 
as degenerate, the sum of the transition amplitudes from these resonances can 
be obtained through the approach in Ref. ||. The transition amplitude for 
the nth harmonic oscillator shell is 



o n = ol + ol 



(23) 



where the amplitudes 0\ and 0\ have the same meaning as the amplitudes 
Aifj and A4fj in Eqs. [11] and [b|, and we have 



1 

+ 3 



o z 



Vfla 



2m, 



[iA ■ (e x k) 



(A x (e x k))] 



1 + 



u 7 
2m n 



a ■ e H -a ■ Ae 



9a 2 m n 




ke • q 



and 



O 2 =0. 



(24) 



(25) 



Generally, the resonances with large quantum number n become important as 
the energy increases. Furthermore, the higher partial wave resonances with 
orbital angular momentum L = n become dominant, which correspond to the 
resonance Gi7(2190) for n = 3 and i/ig(2250) for n = 4 for the rj photo- 
productions. Indeed, only these higher partial wave resonances can be seen 
experimentally, and this is consistent with the quark model predictions. Thus, 
we simply take the masses and decay widths of these high partial wave reso- 
nances as input in Eq. [Ll]. 



3. The Numerical evaluation 

We shall take the same procedure as that in the calculation of the kaon 
photoproductionflj]. In order to take into account of the relativistic effects, 
the Lorentz boost factor is introduced in the CGLN amplitudes 

q) - -— r / i (^ r k, -j-q) (26) 
where i — 1 . . . 4, and ~r (— r-) is a Lorentz boost factor for the initial (final) 

N E N 

state. 
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The parameters in this calculation have standard values in the quark model, 
which the quark mass m q is 0.34 GeV, and a 2 = 0.16 GeV 2 . The masses 
and the decay widths for the S-channel resonances are taken from the recent 



particle data group [ 20 1. In principle, there is only one parameter a v to be 
determined in the numerical evaluation, in which the wavefunctions of the res- 
onances are assumed to have SU(6)<g>0(3) symmetry. However, one should not 
expect that quark model in the 577(6) <g> 0(3) symmetry limit could provide a 
quantitative description of the rj production, since there should be significant 
configuration mixing|B5|]. In particular, the configuration mixing between the 
states iV( 2 P/vf)| 1 and N( a Pm)\ 1 generates a nonzero contribution from the 
resonance 5n(1650), thus affects the r\ photoproduction in the threshold re- 



gion significantly. The evaluations in the potential quark model [I9j show that 
this mixing is indeed very strong. Therefore, we introduce two parameters, 
Cs tl (1535) and Cg u (i65o), to take into account of the configuration mixing ef- 
fects. The contributions from the resonances 5n(1535) and 5n(1650) become 



o Sl 



(1535) 



c Sl 



(1535)O n(2Pm) i- 



(27) 



and 



O 



Su(1650) 



£511(1650) O n{2Pm)1 



(28) 



where On(^p m )\ is given in Table 2. Furthermore, the U-channel contribu- 
tions given in Eq. [H] represent the result in the 577(6) symmetry limit, which 
correspond to Cs u r\s3S) — 1 an d Csn(i650) = 0. There should be an additional 
U-channel contribution for the general coefficients C*5 n , and it is given by 



M 



v 



-M : 



a 



Si i (1535) 



c, 



Si l (1650) 



1 



P f N -k + 5Mlj2 



rj_ qk 



3a 2 



3a 2 



Mr, 



+ 



UJr. 



Ef + M N E l + il7 



N 



UJ ri 



E f N + M ± 



+ 1 



1 + 



N 



2m n 



e,(29) 



where M5 11 ( 1535 ) w Ms^neso) ~ 1-6 GeV in the U-channel. In principle, the 
coefficients Cs u could be obtained from the potential quark models that re- 
produce the baryon spectroscopy. On the other hand, the coefficients Cs u 
obtained from the fitting procedure could provide an important test to the 
wavefunctions in the potential quark models. Moreover, whether the quark 
model could reproduce the large branching ratio for the resonance 5n(1535) 
decaying into r]N channel is still an open question. Therefore, we adopt two 
approaches in the numerical evaluation. First, the wavefunctions of the reso- 
nances are assumed to be in exact symmetry limit, thus only one parameter 
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a v is needed to fit the data. Second, we treat the coefficients, 6*511(1535) and 
Cs n (1650) j and the total decay width of the resonance S , n(1535), 1^(1535) as 
free parameters, and fit them to the differential cross section data. 



The function minimization routine p4] is used to minimize the least square 
function 

X 2 = E [Xt ~ F ' (a 2 , "" ,a " )]2 > (30) 
i (T x l 

where X{ represents the experimental data, ax l corresponds to the error of 
the data, and Yj( CL\ , . . . , CL n ) is the theoretical predictions with parameters 
Oi, • • • , CL n to be fitted. There are about 150 points of differential cross section 
data up to Eiab = 1.45 GeV from the old data setplj, and recent data by 



Homma et a/ ||22j| and by Dytman et a/|| . More recently, the new experimental 
data in the threshold region from the Mainz group has been published ||. 
This set of data differs significantly from the old set data[f2TJ, and recent 
Bates data|| in the threshold region. Therefore, we shall fit the Mainz data 
separately, and the parameters obtained in these fits are summarized in Table 



3. Although there are also few target polarization data|23|], they will not be 



used in our fitting because these data have large errors and do not present any 
systematic behaviour on the target polarization. 

In Fit 1, we assume that the resonances have exact SU (6) (8)0(3) symmetry, 
and the masses and decay widths of resonances come from the recent particle 



data group|2(| . Therefore, there is only one parameters, a v , left to fit the 
data, and we find 

a v = 0.465 (31) 

by fitting it to the combinations of old dataJH], and recent data from 
Bates0. We present the energy dependence of the differential cross sections at 
Qcm = 50°±5° in Fig. 1, at 6 cm = 90°±8° in Fig. 2, and the energy dependence 
of the total cross section is shown in the Fig. 3. Considering that there is 
only one parameter in the calculation, the overall agreement with the data is 
truly remarkable. Our calculation in the symmetry limit also shows that the 
resonance Sii(1535) is less dominant than the data suggest, and the calculated 
total cross section is significantly larger than the data in « 0.9 — 1 GeV 
region. This suggests that the resonance Sii(1650) also plays a significant role 
in addition to the dominant presence of the resonance Sii(1535). Therefore, 
we treat the coupling constant a v , two coefficients Cs n , and the decay width 
rs u (i535) as a free parameter in Fit 2, the resulting parameters are shown in 
Table 3. The resulting fits are also shown in Figs. 1, 2, and 3 respectively. 
It is worth to mention that the smaller total decay width r 5-^(1535) = 0.111 
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GeV is largely due to the recent Bates data|| at E [ab = 0.729 GeV. Since 
only the differential cross section data are used in our fits, the recent total 
cross section data[[7| from ELSA are not used here. Clearly, the threshold 
region is quite crucial in determining the mass and width of the resonance 
5ii(1535). The resonance 5*11(1535) becomes more dominant in this fit, and 
we find a small but negative contribution from the resonance Sii(1650). This 
is qualitatively in agreement with the recent calculation by the RPI group [|i~3| , 
in which the effective Lagrangian approach is used. On the other hand, the 
coupling constant a v is significantly reduced from 0.465 to about 0.14, this 
shows how strongly dependent of the coupling constant a v on the behaviour 
of the resonances Sii(1535) and the resonance Sii(1650). The physical reason 
behind the large reduction of the coupling constant a n is that the threshold 
region is dominated by the resonance Sii(1535), which accounts nearly 90 
percent of the total cross section, thus, a small variation in the Sii(1535) will 
lead to a larger change in the contribution from the Born term. This shows 
that the threshold region alone is not a reliable source to determine the rjNN 
coupling constant a v . 

In Fit 3, we are concentrating on the recent published data from the Mainz 
group 0, in which more systematic differential cross section data are presented 
from Ei a b = 0.716 to 0.788 GeV. It should be pointed out that these data 
are significantly larger than the previous data in threshold region |2T], |5], 22 1, 



therefore, further experimental confirmation is needed. Because this set of data 
is concentrated on the region from the threshold to the mass of the resonance 
5*11 (1535), one could not obtain any reliable information on the contribution 
from the resonance Sii(1650). Thus, we exclude the contribution from the 
resonance S , n(1650) by setting the parameter Cs 1L (ie5o) = 0.0, which is the 
same as that in Fit 1. Thus, three parameters, a v , Cs n (is35) and r Sll (i 53 5), are 
fitted to the data. The calculated total cross section and the data are presented 
in Fig. 4, the agreement with the data is excellent. Moreover, the resulting 
decay width r^n^A is found to be 0.198 GeV, and in very good agreement 
with the simple Breit- Winger fit||. This provides an important consistency 
check of the model. One could also see the possible contribution from the 
resonance S , n(1650) at E iab « 0.8 GeV, where the data suggest that the total 
cross section starts to decrease. More systematic data in Ei ab = 0.8 ~ 1.0 GeV 
region are needed in order to learn more on the structure of both Sii(1535) 
and Sii(1650). 

To highlight the importance of the data in the threshold region in deter- 
mining the properties of the resonance Sii(1535) and coupling constant a v , 
we show the differential cross section at Ei ab = 0.729 GeV in Fig 5, and 
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Eiab — 0.752 in Fig. 6. Notice that there is a significant difference between the 
data from Bates || and Mainz|Q] at E iab = 0.729 GeV, it leads to the change 
from r Sll( i53 5) = 0.111 GeV in Fit 2 to 0.198 GeV in Fit 3. Resolving this 
difference in the future experiments is crucial in understanding the structure 
of the resonance Sii(1535). In the same time, the parameter a n is changed by 
a factor of 3. In fact, the calculation by the RPI group |13[ has shown that one 



could obtain a good fit to the data in the threshold region for a wide range of 
coupling constant a v . More systematic data beyond the threshold region are 
calling for, particularly in the region E^f, = 1.15 to 1.45 GeV, in which there is 
no resonance dominant so that contribution from the Born term becomes rela- 
tively important. It is interesting to note that the calculated differential cross 
section in Fit 2 is in good agreement with the Bate data|| at Eiab = 0.729 
GeV, but smaller at Ei ab = 0.752 GeV, while the results in Fit 3 give excellent 
fits to the Mainz data|| in both cases. 

We should also point out that the resonances Pn(1710) and Pi3(1720) 
also play quite important role in the region Ei ao = 0.9 ~ 1.1 GeV, and one 
should not expect that the quark model in the symmetry limit would provide 
a quantitative description of these resonances. One could also study these 
resonances by inserting coefficients in front of their CGLN amplitudes and 
fitting them to the data. Our calculation provides a framework to study the 
resonance contributions in the meson photoproduction with less parameters. 
This will be investigated in the future with more accurate data in this region. 

The calculation of the target polarizability has also been done with the 
parameters in each fit. We found consistent small polarizabilities at 6 cm = 90° 
from the threshold to Ei ao = 1.0 GeV. If the polarization is indeed large in 
this region as the data suggest [p3| , it might be the evidence that the t-channel 
meson exchange is required. 



4. The Structure of the Resonance 511(1535) 

There were probably two important motivations to study the 77 photoproduc- 
tion in the threshold region; to determine the 7]NN coupling constant a v and 
to study the structure of the resonance Sii(1535). Our calculation shows that 
the threshold region might not be a reliable source to determine the coupling 
constant a eta because of the dominance of the resonance Sii(1535). However, 
one might be able to learn more about the structure of the resonance 5n(1535) 
in the quark model. The study|L3| by the RPI group shows that one could 
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determine the quantity £ from the 77 photoproduction, which is defined as 



£ = JxT v Ai/T T 



(32) 



where x' = M^k/qM^. One can obtain an analytic expression of the quantity 
£ from the CGLN amplitude in Table 2, which is given by 



\ 



a v a e ir(Ef + M N ) C Sll (i 5 35) 



Ml 



6r 9 



2uJ r 



2q2 

3a 1 



UJ r . 



Ef + M N 
Ikl \ 



+ 1 



2m n 



_qf+kf 

e 6a 2 



(33) 



Because q 2 < 1, the quantity £ for the resonance S , n(1535) is not sensitive to 
the parameter a 2 related to the internal structure of the baryon wavef unctions. 
After including the Lorentz boost factors in Eq. EB, we obtain 



0.186 From Fit 1 
0.208 From Fit 2 
0.220 From Fit 3 



(34) 



in the unit of GeV -1 . This is indeed in good agreement with the result 
£ = 0.22 ± 0.02 GeV -1 in ref. |13[. It is not surprising that the quantity 
£ in Fit 1 is smaller, because the resonance Su (1535) is less dominant in 
the symmetry limit than the data suggested. Therefore, assuming that the 
branching ration for the resonance 5n(1535) — > 7]N is around 0.5, we have the 
helicity amplitude; 

81 From Fit 1 

78 From Fit 2 (35) 
t 111 From Fit 3 

in the unit of 10~ 3 GeV - ^, which are very consistent with the results of Ref. 
TBI in Fit 1 and 2, and of Ref. M in Fit 3. 



The advantage of the quark model calculation is that the helicity amplitude 

A\ and the decay width r5 11 (i 5 3 5 ) can be predicted separately. In the symmetry 

2 

limit, the helicity amplitude A\ is given by |26| 



A\ 



2a P 7ruA 



1 

3a 



2m n 



(36) 
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and the decay width in the S , n(1535) 
the coupling constant a^; 



rjN channel is expressed in terms of 



511(1535) 



( V N) 



ar,(E + M N )\q 


a 2 1 


' u v q 2 




2M R 


M%\ 


m q 3a 2 


[E + M N J J 



q 



We have 



and 



A\ = 148 10" 3 GeV^i 



(37) 
(38) 

r Sll( i535)(^) = 23.4MeV (39) 

after including the Lorentz boost factors. Comparing this predictions with 
the results in Fit 1, the helicity amplitude A P i predicted by the quark is twice 

2 

larger than the data, while the decay width is about factor 3 smaller. The fact 
that there is a factor 2 between the old data and the quark model calculations 
for the helicity amplitude A\ has been known for some time[|IO|, [J, and it was 



speculated that this might be an indication of the configuration mixing [pEf . 
However, the systematic calculations with the configuration mixings in the 
Isgur-Karl model shown [26|] that the configuration mixing effects are unable 
to reduce the helicity amplitude A\. Therefore, it is particularly interesting 

2 

that the new data set from Mainz has bring the helicity amplitude A\ in Fit 

2 

3 much closer to the quark model predictions; the simple Breit-Wigner fit in 
Ref. || also gives 

A? = (125 ±25) l(T 3 Ge\H, (40) 

2 

which is even closer to the quark model result in Eq. [38| 

On the other hand, the large r]N branching ratio for the resonance Su (1535) 
has not been fully understood. The calculation by Koniuk and IsgurfllQfl also 
shows that the transition amplitude for 5*11(1535) — > r]N is about 50 percent 
smaller than the data, whose calculation is similar to this approach. This is 
consistent with the fitted coefficient Cs u (1535) ~ 1.5 in both Fit 2 and 3, and it 
is unlikely that this can be explained by the configuration mixing effects. One 
of the effects that has not been taken into account in this study is the finite 
size of T) mesons, which is the one of the major motivations of the quark pair 
creation model|27|]. The calculation has been partly done in Ref. p8| , which 
gives a much larger decay amplitude. The problem is that the r]NN coupling 
that can also be obtained in this approach was not given, thus there is no base 
to judge if the parameter used in the calculation is reasonable. 
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The understanding of the T]N branching ratio of the resonance S f n(1535) 
may be the key to its underlying structure. It has been discussed for some time 
in the literature that the resonance A (1409) might be a Kaon-nucleon binding 
state, whose mass is just below the kaon-nucleon threshold. Notice that the 
mass of the resonance Sii(1535) is just below the threshold of Kaon production, 
7 AT — ■> KA and 7 A" — > A'E, it would be interesting to study the possibility 
that the resonance Sii(1535) is a combination of g 3 and KA or KT, binding 
states, which suggests that the state A(1409) as a KN binding state might 
not be an isolated case. The threshold behaviour of Kaon photoproduction, 
7 A" — > KA and 7 A" — > KH, may provide us further information in this regard, 
because it is dominated by the S wave resonances and the Born terms. 

5. The Conclusion 

The first quark model calculation is presented for the 77 photoproduction, which 
provides a very good description of the rj photoproductions with less parame- 
ters. We show that the threshold region is not a reliable place to determine the 
7]NN coupling constant, which strongly depends on the properties of the res- 
onance Sii(1535). One should extend the study to Ei a b = 1.4 GeV region, in 
which no resonance is dominant, so that the contribution from the Born term 
could be determined more reliably. If there is any indication from the recent 
Mainz data, it might be that the old set of the data may become irrelevant. 
Certainly, the future experiments planned at various facilities, in particular at 
CEBAF, will provide us much more accurate information on the rj photopro- 
duction that will reach beyond the threshold region. The results here show 
that the quark model approach will certainly be a very effective tool of study- 
ing the underlying structure of baryon resonances from the photoproduction 
data. 
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Table 2: The CGLN amplitudes for the S-channel baryons resonances for the 
proton target in the SU(6) <g> 0(3) symmetry limit, where k = |k|, q = |q|, 
and x = ^. The CGLN amplitudes for the iV( 4 P M ), iV( 4 S M ), and N{ A D M ) 
states are zero due to the Moorhouse selection rule, see text. 

States /i / 2 / 3 / 4 



U) 7 
3o^ 



H 2 Pm)\- ^(i + ^J o 

^( 2 Pm)§- + 

o o o 

iv V-^s; 2 36« 2 \ lr 2m,j 216m,a 2 36a 2 U 

/Vl' 2 n'|5+ fcV / -i , fc_\ A: 2 ( 2 lA k k 2 x 

1 "K ±y s) 2 180a 2 I 1 T 2mJ 144m g a 2 V 5^ 180a 2 36ga 4 

^( 2 ^m)| + o ^ 

A^( 2 -Dm)| + tefrfi + db) 216m„a 2 3^ 



/Vr 2 n + fc V (~\ I fc 2 ("„2 lA k k 2 x 

i\yis M ) 2 180a 2 V 2m, y 144m„a 2 V 5 J 180a 2 36ga 4 



Table 1; Meson transition amplitudes A in the simple 
harmonic oscillator basis. 

(N, L) Partial Waves A 

p + 

s ^-(wS^ + i)iS 

(2,0) P Sf-(^ + l)S 

(2.2) P ^-(^ + l)S 

(2,2) F -(i^ + l)S 



20 



Table 3; The parameters obtained from 
different fits. The decay width r 5^(1535) 
has a unit of GeV. 





Fit 1 


Fit 2 


Fit 3 




0.465 


0.139 


0.435 


Csii(1535) 


1.0 


1.510 


1.608 


CSil(1650) 


0.0 


-0.036 


0.0 


rs u (i535) 


0.150 


0.111 


0.198 



Figure Captions 



1. The energy dependence of the differential cross section at 6 cm = 50° ±5°. 
The solid line represents the result from Fit 1, and the dash line from 



Fit 2. The data come from Refs. [ETl E2 



2. The same as Fig. 1 at = 90° ± 8°. 

3. The energy dependence of the total cross section, the solid is the result 
from Fit 1, and dash line from Fit 2. 

4. The result for the total cross section from Fit 3, in which the parameters 
are fitted to the data from Mainz ||. 

5. The differential cross section at E iah = 0.729 GeV. The solid line rep- 
resents the result from Fit 2, and dash line from Fit 2. The data come 
from Ref. (triangle) and Ref. |§ (circle). 

6. The same as Fig 4 at E lah = 0.752 GeV. 
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